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ABSTRACT 



Aims. We present an overview of a high-mass star formation region through the major (sub-)mm, and far-infrared cooling lines 



, to gain insight into the physical conditions and the energy budget of the molecular cloud. 
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Methods. We used the KOSMA 3m telescope to map the core (10' x 14') of the Galactic star forming region DR21/DR21 (OH) 
in the Cygnus X region in the two fine structure lines of atomic carbon (Cl 3 Pi— 3 Po and 3 P2— 3 Pi) and four mid- J transitions 
of CO and 13 CO, and CS J = 7—6. These observations have been combined with FCRAO J = 1—0 observations of 13 CO and 
C ls O. Five positions, including DR21, DR21 (OH), and DR21FIR1, were observed with the ISO/LWS grating spectrometer 
in the [O i] 63 and 145 fim lines, the [C n] 158 /jm line, and four high- J CO lines. We discuss the intensities and line ratios 
at these positions and apply Local Thermal Equilibrium (LTE) and non-LTE analysis methods in order to derive physical 
parameters such as masses, densities and temperatures. The CO line emission has been modeled up to J = 20. 
O Results. From non-LTE modeling of the low- to high-J CO lines we identify two gas components, a cold one at temperatures 
£j of T kin ~ 30 - 40 K, and one with T kin ~ 80 - 150 K at a local clump density of about n(H 2 )~ 10 4 - 10 6 cm" 3 . While the 
^ cold quiescent component is massive containing typically more than 94 % of the mass, the warm, dense, and turbulent gas is 
• . dominated by mid- and high-J CO line emission and its large line widths. The medium must be clumpy with a volume- filling 
^ ■ of a few percent. The CO lines are found to be important for the cooling of the cold molecular gas, e.g. at DR21 (OH). Near 
k> the outflow of the UV-heated source DR21, the gas cooling is dominated by line emission of atomic oxygen and of CO. Atomic 
. and ionised carbon play a minor role. 

_ _ i 

Key words. ISM: clouds - ISM: abundances - Radio lines: ISM - Line: profiles - Stars: formation - ISM: individual objects: 
DR21, DR21(OH) 



1. Introduction 

Emission lines of C n, C I, and CO in molecular clouds 
are among the most important diagnostic probes of star 
formation and particularly important for studying the in- 
fluence of massive stars and their UV radiation on the 
environment. Obtaining the physical parameters of such 
regions is often done with simplified radiative transfer 
assumptions. Local Thermal Equilibrium (LTE) and the 
Large Velocity Gradient (LVG) analysis have proven to 
yield good r esults for the bulk of the cold gas (see, e.g., 
iBladd pOOOj) ). However, the assumption of a homogeneous 
medium does not represent the typical physical conditions 
of a high mass star formation region. Instead, the observ- 
ing beam will contain several gas components with differ- 
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ent excitation conditions and with unknow filling factors. 
The observed line transitions are sensitive to gas at dif- 
ferent density, temperature and opacity, corresponding to 
different regions. By selecting designated lines as tracers 
of the different parts of the gas, one can obtain a rather 
complete coverage of excitation conditions accounting for 
different physical and also chemical conditions within the 
beam. A more differentiated analysis using non-LTE ra- 
diative transfer algorithms (see lvan Zadelhoff et alJ l|2002h 
for an overview) is needed for this purpose. It can provide 
a more detailed description of the structure of the cloud in- 
cludi ng dumpiness and temperature or densi t y gra dients 
(e.g- IOssenkonf et alJ l)200lh. IWilliams et all l)2004l) i. 

Using the results of this analysis, we can try to retrieve 
a self-consistent picture including all major processes driv- 
ing and driven by high mass star formation (e.g, with 
respect to turbulence or chemistry). Thus, we obtain a 
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Fig. 1. View of the DR21/DR21 (OH) high mass star forming region. Left: Overlay of the KOSMA 13 CO 3- 2 map (cf. 
ISchneider etaTl l|2006|) / ISBS200^ in contours on a Spitzer IRAC 8^m greyscale map l|Marston et alJl2004j) . Contours 
are on a T m b scale and range from 16.8 to 100.8 in steps of 16.8 Kkms -1 (a = 5.6Kkms _1 ). The integration interval 
is —10 to 20 kms -1 The left-top circle indic ates the HPBW o f KOS MA. The map also shows 6 positions of radio 
continuu m sources identified as OB stars bv iRoelfsema et all l|l989|) (white crosses), 2 infrared sou rces and maser 
posit ions ijGarden et aljfl986|) (small white circles), and 8 compact 1.3- mm dust emission sources l|Chandler et alJ 
1993) (open rectangles). Right: Overlay of the Spitzer map at 4.5 um with contours of the KOSMA 12 CO 3—2 map 
(integration interval of —20 to 30 kms -1 ; discussed in SBS2006) ranging from 18 to 450 in steps of 48Kkms~ 1 
(er = 6Kkms~ 1 ). The labeled circles mark the positions and beam sizes of the ISO/LWS observations that we use 
to constrain the physical paramters as discussed in the text. The LWS beam is comparable to the resolution of the 
KOSMA maps. 



detailed view of the Photon Dominated Regions (PDRs) 
in the environment of young stars. PDRs are transition 
regions from ionized and atomic to dense molecular gas 
in which physics and chemistry of the gas are dominated 
by the UV radiation from young massive stars. Processes 
heating and cooling the clouds are qualitatively well un- 
derstood on a theoretical level (e.g . , iTielens fc Hollenbachl 
ll985tlSternberg k. Dalgarndll989|) . and in general we find 
a reasonable agreement between the theoretical predic- 
tions and observations of the brightest cooling lines. Still, 
there is an ongoing debate about the details of the en- 
ergy input contributing to the heating. We find a combina- 
tion of mechanisms such as stellar far- ultraviolet radiation 
(FUV), dynamical energy through shocks or supernovae, 
or an enhanced cosmic ray rate. There are also many 
open questions with respect to the chemical composition 
of these regions. The role of nonequilibrium chemistry is 
still not well understood. For instance Polycyclic Aromatic 
Hydrocarbo ns (PAHs) as a r e servoir of fresh gas phase 
carbon (e.g-. IOka et alJ l)2004|) . lHabart et al.l (|2005J)) may 
change the overal picture. 



Emission from warm carbon in the form of atoms, ions 
or in carbonated molecules (mainly CO) provides a sig- 
nificant fraction of the gas cooling at sub-millimeter and 
far-in fr ared (FIR) wave l engths ( e.g..lGiannini et al.l (12 000. 
20011). ISchneider et al] d2003h . iKramer et a.l J ll2004l) . or 
Bradford et all l|2003H2005h ). Finding the prevailing cool- 
ing channels among these and other species (e.g., [Oi], 
H2O, OH) may provide significant insight into the de- 
tails of the chemical composition and the energy balance 
of PDRs. Since systematic observations of sub-mm- and 
FIR-line cooling lines with a good spatial resolution are 
still sparse, data to test the model predictions are still 
largely missing and many questions concerning the details 
of the chemistry and energy balance remain unanswered. 
By studying a well-known, massive Galactic star forming 
region, we aim at gaining new insight into these problems. 

The DR21 region 

The DR21 Hn-region/molecular cloud is part of the 
Cygnus X complex of molecular clouds. Cygnus X 
is known to host a large number of H n-regions 
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Table 1. List of observational parameters for KOSMA. The columns give line frequency, receiver used, main beam 
efficiency 7? m b, half power beam width, observing mode, number of mapped positions, velocity channel width Av (the 
resolution is a factor of ~ 1.4 larger), mean atmospheric zenith opacity (T(f tm ), and observing period. 



Species 


Transition 


Frequency 
[GHz] 


Receiver 


Imb 


HPBW a 
["] 


Mode 


No. 


Av 
[m s" 1 ] 


(r at m) 


Observin 


g period 


13 C q6 


J = 3-2 


330.5880 


DualSIS 


0.68 


80 


OTF c 


693 


308 


0.26 




01/2006 


12 c 32g 


J = 7-6 


342.8829 




0.68 


80 


OTF c 


525 


294 


0.08 




12/2005 


12 CO 


J = 3-2 


345.7960 




0.68 


80 


OTF c 


525 


294 


0.08 




12/2005 


12 CO 


J = 4-3 


461.0408 


SMART 


0.5 


57 


OTF c 


3458 


677 


0.85 - 1.4 




1/2005 


[ 12 ci] 


3 Pi- 3 Po 


492.1607 




0.5 


55 


DBS d 


1380 


630 


0.5 - 1.6 


1/2003 - 


- 3/2004 


13 co 


J = 6-5 


661.067 




0.4 


~40 


PSw c 


33 


307 


1.0 - 1.3 




1/1998 


12 co 


J = 6-5 


691.473 




0.4 


~40 


PSw c 


58 


295 


1.0 - 1.3 




1/1998 


12 co 


J = 7-6 


806.6518 




0.31 


42 


DBS d 


964 


775 


0.4-1.6 


1/2003 - 


- 3/2004 


[ 12 ci] 


3 P 2 - 3 Pi 


809.3420 




0.31 


42 


DBS d 


964 


775 


0.4-1.6 


1/2003 - 


- 3/2004 



a Th e main beam efficiency ?7 m b and the HPBW were determined from cross scans of Jupiter. 

b See ISchneider et alJ (|2fl06) for a 13 CO 3—2 map covering a larger part of Cygnus X. 

c OTF - On-The-Fly mapping mode with periodic switching to a emission free reference position at 

q(J2000) = 20 h 37 m 10 s <5(J2000) = 42°30'00". 
d DBS - Dual Beam Switch observing mode with a chopping secondar y mirrow in a zimuth. 
e The CO and 13 CO 6-5 observations were carried out on a 35" grid <IKosteilll998t) in a Position-Switched 

observing mode (PSw). The OFF position was at a( J2000)=20 h 26 m 43 8 .56, S( J2000)=42°09'39".5. 



l|Downes fc Rinehartll96rilWendkerll984 associated with 
mol ecular clouds as seen e.g., in the CfA 12 CO 1—0 sur- 
vey ijLeung fc Thaddeusl l|l992|) '). An extended (11 square 
degrees) map of the Cygnus X region in the 13 CO 2—1 
line and smaller maps in the CO and 13 CO 3—2 lines were 
obtained with the KOS MA telescope and are pr esented in 
ISchneider et"ai] l)2006|) (from now on lSBS2006j) . The dis- 
tance to DR21 was sometimes estimated to be ~ 3kpc, 
base d on visual interstellar extinction and radial velocity 
(e.g. , ICampbell et all 1982). For this paper, we follow the 
argumentation from SBS2006 based on a morphological 
comparison of the molecular line data with mid-IR emis- 
sion and the OB associations, and favour a value of 1.7 kpc 
as distance toward DR21 and DR21 (OH). However, for 
easier comparison with the literature, we explicitly give 
distance depended physical properties as a function of dis- 
tance. 

The DR21 region has been subject to numerous stud- 
ies at different wavelengths. Several active star forma- 
tion sites and cometary shaped H n-regions w ere identified 
in maps of radio continuum e mission (e.g., iDickel et alJ 
1983[ |Cyganowski et alj| 2003tl . Ammonia observations by 
Mauersbereer et al lll985h were interpreted as evidence 
for high densities (~10 5 5 cm~ 3 ) paired with a high de- 
gree of dumpiness. H2O masers l)Cenzel fc Downeslll977l) 
and a recently obtained map of 1.3- mm continuum emis- 
sion l|Motte et alJ l2005|) show that DR21 belongs to a 
north-south orientated chain of massive star forming com- 
plexes. Additional indication of active star formation is 
also given by observations of vibrationally excited H2, 
tracing the hot, s hocked gas in t he vic inity of DR21 and 
DR21 (OH ) fe.g.. iGarden et all jl986|) . and more recent 
results by iDavis et al.l l|2006f) ). H 2 images reveal an ex- 



tended (5'=2.5pc(D/1.7kpc)) north-east to south-west 
orientated outflow centered on DR21. To the west, the 
outflow is highly collimated, whereas toward the eas t 
side, gas is expan ding in a blister style ijLane et alJ fl99C0. 
Ijaffe et all (^L989) found very intense CO 7—6 emission 
peaking on the H 11-region and following the H2-emission. 

In this paper, we present the results of sub-mm map- 
ping observations with KOSMA and discuss far-IR line in- 
tensities observed with ISO. This comprehensive data set 
allows to study the spatial and partly also the kinematic 
structure of all major cooling lines of the ISM in the vicin- 
ity of DR21 and DR21(OH). In Sect. 3, we analyse the 
line emission and continuum data toward five positions, 
by comparing the spatial distribution, line profiles, and 
the physical properties through selected line ratios and 
the far-IR dust continuum. In Sect. 4, we model the in- 
tegrated line intensity ratios and absolute intensities with 
radiative transfer models. The impact of the line emis- 
sion on the total cloud cooling is dicussed in Sect. 5. A 
summary is given in Sect. 6. 

2. Observations 

2.1. KOSMA 

The region centered on DR21/DR21 (OH) was mapped 
in the two atomic carbon fine structure transitions at 
492GHz (609 /an, 3 Pi- 3 P ; hereafter 1-0) and 809 GHz 
(370 jtrai, 3 P 2 - 3 Pi; hereafter 2-1) and the mid-J ro- 
tational transitions of CO (J=3— 2, J=4— 3, J=6— 5, 
and J=7-6), 13 CO (J=3-2, and J=6-5) using the 
dual channel 345 & 230/660 GHz SIS receiver, and 
the SubMillimeter Array Receiver for Two frequencies 
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Fig. 2. Integrated intensity map of CO 4-3 (left) and CO 7-6 (right) of DR 21 & DR21 (OH). The integration interval 
is —20 to 30 kms^ 1 . The contours for CO 4—3 run from 66(3cr) to 462Kkms~ 1 in steps of 3(7, and from 68(2<r) to 
671 Kkms- 1 in steps of 2a for CO 7-6. The HPBW is indicated in the upper left corner. The positions of ISO/LWS 
observations (cf. Fig. ^) are indicated by circles. The dashed white lines correspond to the cuts shown in Fig. 



Table 2. Positions of ISO observations in absolute coor- 
dinates. DR21 C is centered on the Hn-region. Projected 
distance of DR21 E (resp. DR21 W) from position DR21 C 
is 0.77 pc (resp. 0.89 pc) at a distance of 1.7 kpc. The TDT- 
number is given in the second column. 



Name 


TDT-No. 


R.A. [J2000] 


Dec. [J2000] 


DR21 E/East 


15200785 


20 h 39 m 09M4 


42°20'05".0 


DR21 C/Center 


15200786 


20 h 39 m 00 s .93 


42°19'42".0 


DR21 W/West 


15200787 


20 h 38 m 52 s .23 


42°19'00".5 


DR21 (OH) 


34700439 


20 h 39 m 00 s .95 


42°22'37".6 


DR21 FIR1 


35500317 


20 h 39 m 00 s .2 


42°24'34".6 



(SMART) at the KOSMA-3m telescope. SMART is a 
dual-frequency, eight-pixel SIS-heterodyne receiver that 
simultaneously o bserves in the 6 50 and 350 (im atmo- 
spheric windows ijGraf et al.ll2002() . All observations were 
taken in the period between 2003 and 2006. 

An area of ~ 10' x 14' was mapped in the two [C i] lines 
and the CO 7-6 line on a 27.5" x 27.5" grid with a total 
integration time of 160 s per position. The observations 
before 2005 were performed in Dual-Beam-Switch mode 
(DBS), with a secondary mirror chop throw of 6' fixed 
in azimuth. Due to the extent of the DR21 cloud com- 
plex, self-chopping is present in some of the DBS-mode 
observations. In Appendix [5] we describe how we pro- 
cessed these data in order to reduce the impact of the 
self-chopping. CO and 13 CO 3-2, and CO 4-3 were ob- 



served in an On-The-Fly (OTF) mapping mode on a fully- 
sampled grid. A zeroth order baseline was subtracted from 
the calibrated spectra. In addition to CO 3—2, we simulta- 
neously obtained the CS 7—6 line in the image sideband of 
the 345 GHz channel. Pointing was frequently checked on 
Jupiter, the sun, and then confirmed on DR21 itself. The 
resulting accuracy was better tha n ±20". L arger 13 CO 
3—2 and 2—1 maps are presented in SBS2006. Due to their 
lower spatial resolution, we did not consider the J=2— 1 
data here. 

The 12 CO an d 13 CO 6—5 observations were presented 
in iKosterl l|l998h . The typical DSB receiver temperatures 
in the high frequency channel was between 160-220 K. The 
13 CO 6—5 observations cover a region of ~ 3' x 3' centered 
on DR21 while the 12 CO observations extend further to 
the north covering a region of ~ 3' x 6'. 

Atmospheric calibration was done by measuring the at- 
mospheric emission at the reference position to derive the 
opacity. Sideband imbalan ces were corrected using stan- 
dard atmospheric models l|Cernicharolll985|) assuming a 
receiver sideband gain of 0.5. Spectra taken simultane- 
ously in the two frequency bands near 492 and 810 GHz 
were calibrated with an averaged water vapor value. The 
[Ci] 2—1 and CO 7—6 lines, observed simultaneously in 
two sidebands, were also corrected for sideband imbalance. 
Their pointing is identical and the relative calibration er- 
ror between these two transitions is very small. We scaled 
the antenna temperature data to the T m b scale using main 
beam efficiencies r] s as listed in Table ^ The HPBWs 
(40" to 80") and r^mb were determined using continuum 
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Fig. 3. Integrated intensity map of [Ci] 1—0 (left) and of [Ci] 2—1 (right). The integration interval is —10 to 20 
kms -1 . The contours for Ci 1—0 run from 12(2er) to 60Kkms~ 1 in steps of 2a, and from 20(1ct) to 108Kkms _1 in 
steps of a for [Ci] 2-1. 



scans on Jupiter. From observations of standard calibra- 
tion sources, we estimate that the absolute calibration is 
accurate to within ~ 15%. Most spectra were taken under 
excellent weather conditions with typical values of the at- 
mospheric opacity along the line-of-sight better than 1.0 
at 350 and 600 //m, and better than 0.5 above 870 /im. 



2.2. ISO LWS Observations and Data Reduction 



43°30'00" 



ISO L ong Wavelength Spectrometer (LWS, IClegg et"aT] 
l|l996|) 1 43 - 197 /im grating scans (AOT L01) were ob- 
tained for 5 positions in the DR21 region from the ISO 
Data Archive (IDA) (see Tableland Fig.0). We removed 
duplicate scans, overlapping positions, or data with a very 
poor spectral resolution from the sample. 

The grating scans contain the FIR continuum, the 
atomic fine structure lines of [O i] at 63 and 145 fim, 
ionised atomic lines of [Cll] at 158 /im, [Olll] at 52 and 
88 fim, [Nil] at 122 /im, and the high- J CO molecular lines 
J=14— 13 through J=17— 16. All spectra were processed 
with the ISO Spectral Analysis Package (ISAP, v. 2.1). In 
ISAP, the data were deglitchcd by hand, defringed (detec- 
tors 4 — 9), and corrected for flux clipping (the extended 
source correction). Some artefacts remain after the pro- 
cessing (see Fig. and were not corrected. Line strengths 
were measured within ISAP. See Table 01 for further de- 
tails. 
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Fig. 4. Integrated intensity map of CS 7—6. The inte- 
gration interval is —10 to 4 kms^ 1 . Contours run from 
1.8(3cr) to 13.0 Kkms- 1 (in steps of 3a). 



2.3. FCRAO 

We use 13 CO and C 18 1—0 data from a large scale map- 
ping project of the whole Cygnus X region using the 
FCRAO 14m telescope. These data were obtained be- 
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Fig. 5. Spectra of all observed species at the five ISO positions. Grey filled spectra outlined in black: observed lines, 
thin grey lines: line profiles modeled with the radiative transfer code SimLine (see Sect. l4.2|) . All spectra are shown at 
a common spatial resolution of 80" and on a T m b scale. The baselines are shifted in steps of 30 K. [Ci] was multiplied 
by a factor of 2, CS 7—6 and C 18 1—0 by a factor of 4. J = 1—0 transition data are FCRAO observations. The 
dashed lines indicate the central velocity of —3 kms" 1 , and the dotted lines the secondary velocity component at 
about 9 kms -1 , which was excluded from the fitting. 



tween 2003 December and 2006 January and will be pre- 
sented in more detail in upcoming papers (Simon et al., 
in prep.). The data were observed using the 32 pixel array 
SEQUOIA in an On-The-Fly mapping mode (at an angu- 
lar resolution of —50"). In this paper, we use smoothed 
data to have the same angular resolution of 80" as the 
KOSMA data. 



8 /im emission therefore highlights the photon-dominated 
surface regions of molecular clouds (PDRs). The molec- 
ular outflow along DR21 E, C, and W is well traced by 
4.5 fim emission, mainly due to strong lines of shock ex- 
cited or UV-pumped H2 v=0— S(9). No ro- vibrational 
lines of CO v =1— and only w eak Bra have been found 
in this band {Smith et alJl2006|) . 



3. Results 

Figure n gives an overview centered on DR21 (OH) (des- 
ignated as W7 5 S by some authors) . The Spitzer 8 ^m and 
4.5 /zm data bv lMarston et all l|2004|) show the small scale 
structure of this region at angular resolutions between 
1.3"~3". Several bright IR objects, in particular near 
DR21 C 1 , delineate the dense ridge of gas and dust in both 
wavelength bands. The map of 8 /im emission is dominated 
by many streamers of emission, which appear to connect 
to the ridge, and diffuse emission. This is mainly emission 
from Polycyclic Aromatic Hydrocarbons (PAHs), which 
are the prime absorbers of far-ultraviolet (FUV) radiation 
from the surrounding and embedded massive stars. Bright 

We will use the ISO/LWS position names troughout the 
paper (cf. Table 0. 



3.1. Line-integrated maps 
3.1.1. CO and 13 CO maps 

Figure ^ (contours on right panel) and Fig. [21 show the 
KOSMA maps of integrated emission of the rotational 
transitions CO 3—2, 4—3, and 7—6 which preferentially 
trace the warm and dense molecular gas around DR21 and 
along the north-south orientated molecular ridge. Apart 
from the ridge, the most obvious feature seen in the maps 
is the line emission along the outflow axis oriented NE- 
SW. The main axis of the strong 3—2, 4—3 and 7—6 
emission (FWHM -230", resp. ~ 1.9pc (D/1.7kpc)) is 
tilted by 20° with respect to the right ascension axis 
and correlates nice ly with extended H2 emission lobes 
I Garden et aHllflSfift . The red-shifted CO line- wing emis- 
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Table 3. Line integrated intensities in [Kkms' 1 ] (upper table) and in [1CP 7 ergs -1 cm~ 2 sr _1 ] (lower table) of 
the low- and mid-J CO lines observed (including the 13 CO and C 18 isotopomers) with KOSMA and FCRAO, and 
CS 7—6. In order to prevent the +9 kms -1 velocity component from adding to the integrated intensity, the listed 
intensities were integrated only over the —3 kms -1 component. As the 12 CO lines merge with the second velocity 
component, these spectra were fitted using a 2-component gaussian with the +9 kms -1 velocity component masked 
out. All line intensities are reported for 80" angular resolution (higher angular resolution data have been smoothed). 



Position 


CO 


CO 


CO 


CO 


13 CO 


13 CO 


13 CO 


c 18 o 


[Ci] 


[Ci] 


CS 


Units in [Kkms -1 ] 


3-2 


4-3 


6-5 


7-6 


1-0 


3-2 


6-5 


1-0 


1-0 


2-1 
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DR21 E 


348.7 


283.4 


247.2 
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22.38 


69.4 


43.1 


1.66 


21.9 


28.2 


3.0 


DR21 C 


440.0 


462.6 


434.4 
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70.48 


112.2 


80.9 


8.26 


46.6 


58.4 


13.8 


DR21 W 


230.5 


221.9 


177.7 


223.1 


41.02 


47.6 


33.9 


3.67 


26.9 


28.5 


2.9 


DR21 (OH) 


283.2 


197.4 


162.0 


153.2 


67.92 


91.5 


24.9 


9.03 


39.5 


30.2 


13.6 


DR21 FIR1 


223.8 


193.2 


130.2 


107.3 


59.49 


61.9 


26.0 


9.16 


33.4 


33.0 


3.9 
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CO 


CO 


CO 


CO 


13 CO 


13 CO 


13 CO 


c 18 o 


[Ci] 


[Ci] 


CS 


in [10 _7 ergs _1 cm -2 sr~ 


x ] 3-2 


4-3 


6-5 


7-6 


1-0 


3-2 


6-5 


1-0 


1-0 


2-1 


7-6 


DR21 E 


147.8 


284.6 


837.6 


1429.1 


0.307 


25.7 
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0.023 


26.8 


153.4 


1.3 


DR21 C a 
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464.6 


1471.9 


2568.6 


0.967 


41.5 


239.8 


0.112 


56.9 


317.3 


5.7 


DR21 W 


97.7 


222.9 


602.1 


1200.2 


0.563 


17.6 


100.3 


0.050 


32.5 


154.7 


1.2 


DR21 (OH) 


120.0 


198.3 


548.9 


824.1 


0.932 


33.9 


73.8 


0.122 


48.3 


164.7 


5.6 


DR21 FIR1 


94.8 


194.0 


441.2 


577.0 


0.816 


22.9 


77.0 


0.124 


40.8 


179.5 
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Fig. 6. Cut along the DR21 outflow axis (top) and along 
the ridge (bottom) centered on DR21 C in CO 3-2, 4-3, 
and 7—6 integrated emission. The 80" resolution is indi- 
cated as horizontal line. The dashed vertical lines indicate 
the positions of the ISO/LWS observations. 



sion is stronger toward the west side, whereas the blue- 
shifted emission is more pronounced to the east indicating 
that this side is closer to the observer. 

Morphologically, the CO 3—2 map shows a very sim- 
ilar shape in comparison to the CO 4—3 data. Detailed 
differences can be better seen if the maps are smoothed to 
the same resolution. In Fig. we present two cuts along 
the outflow and along the main ridge from these data. At 
DR21, all CO line integrated intensities are very similar. 
Toward position DR21 E, the emission of CO 3—2 and 4—3 



stays at a high level while the 7—6 emission drops. Toward 
the west, all three line intensities decline in a similar fash- 
ion. Away from DR21, along the N-S cut, the CO 3-2 
emission is dominant, whereas CO 4—3 and CO 7—6 are 
sequentially weaker. Strong CO 7—6 emission is limited to 
the DR21 core. DR21 (OH) is located 180" north of DR21 
and clearly separated from the DR21 region by a gap in 
the ridge. DR21 (OH) is bright in the lower J transitions, 
but appears weaker in the J=7— 6 transition. 



13 CO 3—2 (contours on left panel of Fig. ^) was se- 
lected to map the column density in the region. This 
transition traces the molecular ridge showing two maxima 
of emission at DR21 C and DR21(OH) The tota l mass 
of the ridge is 3.5x 10 4 M Q (D/1.7 kpc) 2 (SBS2006). 13 CO 
weakly follows the outflow from DR21E to W. Following 
the molecular ridge further to the north, an elongated 
emission band roughly coincides in position with three far- 
infrared sources. The sout hernmost sourc e, DR21FIR1, 
was observed by ISO /LWS. Ikarston et all l|2004l) find an 
extremly red object (ERO 3) here, illuminating a shell-like 
halo in the 4.5 /im band. The other two, FIR2 and FIR3, 
coincide with MSX point sources. As the FIR sources, this 
northern region is slightly tilted to the east, indicating flux 
contribution to mainly low- J CO emission from all three 
embedded far-IR sources. 13 CO reaches up to FIR3 which 
marks the edge of the dense ridge to the north. 
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3.1.2. [Ci] maps 

C 1 1—0 and 2—1 appear in good agreement with the shape 
of the molecular ridge as traced by 13 CO. Even though the 
emission around DR21 in [Ci] (see Fig. 01 appears more 
centrally peaked than in the CO maps, the main axis be- 
tween positions E and W is clearly visible. The molec- 
ular ridge is seen to extend from the south east through 
DR21 and through DR21 (OH) 3' northwards. DR21 (OH) 
is visible as elongated, less compact region linking-up to 
the far-infrared sources further north. The DR21 (OH) re- 
gion appears to harbour multiple sources as there is no 
central emission peak. The [C i] 1—0 map reveals a south- 
ern peak and extended weaker emission in the north. This 
substructure may n ot be significant becau se the variation 
is only about 2a. IChandler et alJ l)l993[) identified five 
compact sources in the DR21 (OH) region (indicated in 
Fig. ^) using 1.3- mm thermal dust emission at 11" reso- 
lution. Ci 2—1 is, compared to 1—0, relatively weak to- 
ward DR21 (OH) and suggests a lower temperature than 
for DR21C (see below). North of DR21 (OH), the car- 
bon emission stretches over the three far-IR spots FIR1, 
FIR2, and FIR3. Here, the upper atomic carbon line is 
in good agreement with the overall shape of optically 
thin low- J CO isotopomers, i. e., the C 18 2—1 map by 
IWilson fc Mauersbereerl [jigged , and to some extent our 
13 CO 3—2 map. A more distinct morphology analysis is 
not possible due to the remaining noise and some self- 
chopping residuals. 

3.1.3. CS 7—6 map 

The map in Fig. 0] reveals two condensations at DR21 
and around DR21 (OH). The southern peak is centered 
on the H Il-region and follows the direction of the outflow 
to the west. The FWHM is 133 "x 102". The northern 
complex (FWHM 108 "x 156") is extended and may con- 
tain several unresolved cores (see above). The emission 
extends up to the FIR1 position, but does not cover all 
three FIR spots. The lines of CS 7-6 and C 18 1-0 are 
the only lines in our sample that show brighter emission 
toward DR21 (OH) than toward the southern molecular 
cloud /H H-region DR21 . For CS, this has also been ob- 
served bv lShirlevet~ai1ll2003h . 

3.2. Line profiles 

In the following, we concentrate on the 5 ISO/LWS posi- 
tions (as indicated in Fig. QJ. In Fig. O we present spec- 
tra of 12 CO, 13 CO, and C 18 towards these positions to- 
gether with spectra of Ci, and CS 7—6 between —30 to 
30 kins" 1 . We present all observations at a common res- 
olution of 80". In Table we give the corresponding line 
fluxes. 

The 12 CO emission from the DR21/Hn-region shows 
broad wing emission due to the outflowing molecular gas 
(with line widths of more than 15 kms -1 and peak tem- 
peratures well above 20 K). A self-absorption dip at the 



velocity of —3 kms -1 is most pronounced in mid- J CO 
lines at DR21 C where the peak intensity drops by a factor 
of two at the line center. Toward DR21 (OH), CO shows 
blue-skewed, double-peaked lines with self-absorption. In 
both cases, the velocity of the line reversal nicely matches 
the corresponding emission in 13 CO and Ci, and is likely 
caused by the same cold material in the dust lane. 

With widths of 4 to 6 kms" 1 , the line shapes of 13 CO 
are similar to the [Ci] profiles. The 13 CO 6—5 emission 
(E u /kB — 79.3 K) at all five positions indicates the pres- 
ence of warm and dense gas. At DR21 C, 13 CO 6—5 shows 
broad symmetric wing emission that is underlying a nar- 
row component. The velocity extent of the wings is sim- 
ilar to that in mid-J CO lines. A better signal-to-noise 
ratio would allow a more accurate estimate. Only weak 
line wings are visible in the lower 13 CO 3—2 line since it 
traces primarily the quiescent gas. 

DR21 is associated with molecular gas centered at a 
velocity of —3 kms" 1 . In addition to this, weaker emis- 
sion between ~ 7 and 14 kms -1 is associated with the 
W75N complex approximately 11 pc (D/1.7kpc) t o the 
north-west. Some authors (e.g.. iDickel et ai1ll978t) con- 
clude that both clouds are actually interacting~[sBS2006 
shows channel maps revealing a link between both clouds. 
Diffuse emission in the 13 CO J = 1— and both [Ci] tran- 
sitions, and absorption in the red wing of the low- J CO 
lines (up to J = 4—3) suggest that this material is located 
in front of the DR21 complex. We do not fu rther discuss 
this contribution here and refer to lSBS2006l 

3.3. ISO archive data 

Figure shows the ISO/LWS continuum fluxes together 
with spectrally unresolved lines of [C n] , [O i] , [O in] , 
[Nil], and rotational high- J CO transitons. The following 
Sections discuss the dust and line emission at the individ- 
ual positions. 

3.3.1. Grey body fit to the continuum 

In all five cases, the far-infrared dust-continuum flux - 
with strong emission lines masked out - was x 2 -fitted with 
a single isothermal grey-body model F(Tdust, T dust ) A = 
ClB x (T dust )(l - e ~ Td ^) with fixed solid angle (fi = 1.5 • 
10" 7 sr, assuming that the emission fills the 80" aper- 
ture) and a dust opacity Td us t ~ t v A"^ with = 1.5. 
Spe ctral indices are typicall y found between 1.0 and 2.0 
(see IColdsmith et alJ l(l997l) for a discussion of observa- 
tions and grain models). Although the spectral index is 
known to vary with the cloud type and optical depth, we 
assume a fixed for all 5 positions (but see discussion in 
Sect. |4. 231 . A fixed of 1.5 is compatible wit h grain mod- 
els of IPreibisch et all lll993h. IPollack et all l|l994h . and 
observations, e.g., bv I Walker et alJ ((199(T ). The total in- 
frared dust continuum flux and bolometric luminosity per 
position is derived by integrating the fit-model F\. The 
integration interval (1 — 1000 /urn) covers the expected in- 
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Table 4. The ISO/LWS high- J CO and atomic FIR lines at 80" resolution. The transformation between flux densities 
and surface brightnesses was done using the beam sizes given in Table 5.10 of the ISO-LWS handbook v. 2.1. The 
beam solid angle varies between 8 x 10~ 8 sr and 1.4 x 10~ 7 sr. Because the lines are unresolved, the line widths were 
held to the instrumental line width: 0.29 or 0.60 /im for detectors — 4 and 5 — 9, respectively. Correction factors 
between 0.37 to 0.85 for extended sources have been applied. We assume a 30% absolute photometry uncertainty for 
LWS. For DR21 (OH), numbers in parentheses are more uncertain (S/N between 1 — 3) due to a poorer frequency 
resolution. At DR21FIR1, no detection of CO was possible, so we give 3<r upper limits. [Nil] was clearly found only 
toward the outflow positions E and W. At the other positions we give 3cr upper limits. 



Position 


CO 


CO 


CO 


CO 


[Cn] 


[Oi] 


[Oi] 


[Om] 


[Om] 


[Nil] 


in [10~ 4 ergs -1 cm~ 2 sr _1 ] 


14-13 


15-14 


16-15 


17-16 


158 tzm 


63 fim 


145 fim 


52 /im 


88 lira. 


122 ^m 


DR21 E 


0.25 


0.72 


0.51 


0.33 


8.2 


48.2 


2.9 


3.53 


3.08 


0.435 


DR21 C 


1.75 


2.46 


0.99 


0.63 


8.1 


49.4 


6.15 


11.5 


4.64 


< 0.75 


DR21 W 


0.44 


0.69 


0.42 


0.29 


4.7 


10.1 


1.136 


1.78 


2.52 


0.306 


DR21 (OH) 


(0.45) 


(0.99) 


(0.40) 


(0.50) 


2.3 


2.94 


0.845 


(1.27) 


1.56 


< 2.7 


DR21 FIR1 


< 0.28 


< 0.3 


< 0.28 


< 0.25 


3.6 


1.78 


0.45 


0.54 


1.66 


< 0.6 



Table 5. Results of dust continuum single-component fit. (3 was fixed at 1.5, and SI is 1.5 • 10~ 7 sr. The extinction and 
the opacity are linked through A v = 2.5r v log 10 e ss 1.085r v . The total infrared flux is calculated from Im = J F\d\. 
The luminosity L and mass M are given for 1.7 kpc distance. 



Position 


T v 


A v (ISO) 


N(H 2 ) 
[10 22 cm" 2 ] 


Mass 
[M (D/1.7kpc) 2 ] 


Tdust 

[K] 


Im. 

[ergs -1 cm -2 sr _1 ] 


L 

[1O 3 L (D/I.7kpc) 2 ] 


DR21 E 


20 


21.2 ± 1.9 


2.1 


148 


41.6 ± 0.6 


0.55 ± 0.05 


7 


DR21 C 


112 


121 ± 11 


12.1 


850 


42.5 ± 0.7 


3.4 ±0.3 


45 


DR21 W 


34 


36.9 ±3.5 


3.69 


259 


30.7 ±0.4 


0.18 ±0.02 


2 


DR21 (OH) 


238 


259 ± 43 


25.8 


1814 


31.2 ±0.8 


1.3 ±0.2 


17 


DR21 FIR1 


45 


48.4 ± 4.3 


4.84 


340 


33.3 ±0.5 


0.37 ±0.03 


5 



tensity maxima bet ween 50-80 um and c orresponds to the 
range used, e.g., bv lCalzetti fittil foMHjl . Errors of the fit 
parameters were estimated as described in the following: 
we assume (i) the optical depths Tdust and total fluxes are 
accurate to about ~ 10%, mainly due to the remaining 
calibration artefacts within the LW1— 5 bands, (ii) the 
temperature fit error is less than 2 % but the absolute 
error also depends on the quality of SW1 detector data 
(leftmost data points in Figure 01. 

We find extinctions between 20 m <A„<260 m and tem- 
peratures between 31— 43 K (cf. Table EJ). In contrast to 
the extinction, the variation in dust temperature is small. 
Considering that we use a single-temperature model, these 
temperatures are probably biased by the warmer dust. 
Evidence of cold dust peaking above 100 /jm is visible at 
DR21W, E, (OH) and FIR1 in Fig. □ but disentangling 
this component would require longer wavelength data and 
the contribution to the integrated grey-body flux is rel- 
atively small. The co ld du st is studied in det ail by, e.g., 
IChandler etaTl (|l993[) and lMotte et alJ l|2005[) . We derive 
the molecular hydrogen column density fro m the conver- 
sion N(H 2 )/A„ w 0.94 • 10 21 cm -2 mag -1 l(Bohlin et all 
119781 cf. Table El) and if placed at the distance of 1.7 kpc, 
the molecular gas masses vary between 150 — 1800 M . 

The peak of the beam-averaged dust temperature is 
-43 K at DR21 C. With A v > 100 m , the strongest contin- 
uum sources are DR21 C and DR21 (OH). Due to the high 



temperature towards DR21 C, the flux density is higher 
at short wavelengths, but drops below the flux density 
of DR21(OH) above -200 /zm. For DR21 (OH), we es- 
timate a higher mass t han for the DR 2 1 E-W complex 
(1814, resp. 1257M ). IChandler et"afl l|l993|) assumed 
temperatures between 25 K and 40 K and derived a slightly 
lower total mass of 1650MQ(D/1.7kpc) 2 for the sources 
at DR21 (OH). 



3.3.2. ISO/LWS far-IR lines: Carbon Monoxide and 
Ionised Carbon 

The LWS band covers many important far-IR cooling lines 
and tracers of the photon dominated regime of molecular 
clouds. TableH gives the integrated line intensities of [C n] 
(158 /zm), high- J CO, and [Oi] (63 and 145 fan) derived 
from the ISO observations quoted above. In this para- 
graph, we focus on the carbon bearing species needed for 
Sect. 0] and postpone the discussion of the other lines to 
the end of Sect. El 

At DR21, highly excited CO emission is visible at all 
three positions, but most prominently around the H n- 
region at DR21 C The center position with more than 
2.4 • 10 _4 ergs" 1 cm~ 2 sr _1 in mid- and high- J CO is also 
the brightest continuum source in the mapped region. The 
emission stems from large columns of warm quiescent gas 
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Fig. 7. ISO/LWS Spectral Energy Distribution (SED) of 
the DR21 region (upper panel), and the DR21 (OH) and 
DR21FIR1 region (lower panel). Curves indicate the best 
greybody fit corresponding to the values given in Table 03 



that is heated from an embedded object. Toward the East 
and West positions, where vibrationally excite d H 2 out- 
lines the molecular outflow ijGarden et alj fl986). the FIR 
continuum is somewhat weaker and more importantly, 
the continuum flux does not correlate with the unusually 
strong FIR line emission. This gas is visible in the line 
wings of mid- J CO lines (Fig. |5j) and in bri ght H? line 
emissi on at 4.5 /xm in the IRAC map (Fig. fTf) . lLane et "all 
( 199(J conclude that neither a single C- nor a single J- 
shock can account for the observed far-IR line emission 
alone and additional heating through FUV radiation is 
required. Although DR21 (OH) has the second strongest 
continuum, the sources associated with this position arc 
surprisingly inactive in far-IR line emission. We find four 
high-J CO line detections, but due to a lower scanning- 
resolution, the values given in Tabled have higher uncer- 
tainty. The far-IR source FIR1 is the faintest with respect 
to the line fluxes, although the continuum is stronger than 
towards DR21W. High-J CO was not detected at FIR1, 
so we give upper limits in the table. 

[Cn] is detected at all five positions. The emission 
peaks at DR21 E and toward the West, the flux de- 
creases by almost 50%. The low [Cn] flux of 2.3 • 
10 _4 ergs _1 cm~ 2 sr _1 at DR21 (OH) may indicate a lack 
of ionised carbon or that the temperature is too low for 
excitation in the upper energy level (E u /fcs = 91 K above 
ground). At DR21FIR1, [Cn] is again slightly enhanced 
compared to DR21 (OH). 



4. The physical structure of DR 21 

In this section we analyse the line properties of the ob- 
served carbon species from ISO/LWS data and combine 
these results with the KOSMA observations. We start with 
an LTE approach and then model the line emission with 
full radiative transfer. 



4.1. Comparison of CO, C, and C + 

We estimate the composition of the molecular gas with 
respect to the carbon abundance in the main gas-phase 
species CO, C, and C + . This is done by assuming LTE 
and optically thin line emission. We therefore use the 
less abundant isotopomeres 13 CO and C 18 0. We list the 
results for excitation temperature, column density and 
abundance in Table 



4.1.1. Atomic Carbon 

The C I excitation temperature needed to produce the ob- 
served line ratios can be estimated since both C I lines are 
likely to be optically thin (cf. line profiles in Fig. . The 
line temperature ratio i?ci falls between 0.8 at DR21 (OH) 
and 1.3 around and east of DR21. The temperature is de- 
termined from this ratio using T cx = 38.8K/ln(2.11/i?c : ) 
(Eq. (ESJ). The range (38 K at (OH) to 79 K at E) indi- 
cates the origin of the atom i c carb on lines in a warm envi- 
ronment. [Zm^i^zin^^^^ l|l988h derived a temperature 
of 32 K from the two C I lines at DR21 (OH) which is close 
to our value at this position. Our results for T ex are lower 
compared to other Galactic high mass star forming re- 
gions (e.g., W3 Main and S 106 with typical ly T ny > 110 K, 
iKramer et al"1r2004t Isdmeider et al.ll2003l) . but substan- 
tially higher compared to the a yerage temperature in the 
inner Galactic disk of - 20 K l)Fixsen et alJll999l) . The 
beam averaged column density, 3.2 — 6.8 • 10 17 cm~ 2 (Eq. 



2J|), shows variations by a factor of 2 or less, peaking 
at the position DR21C and around (OH), i.e., towards 
the densest condensations along the ridge. This is remark- 
able since the H 2 column density derived from the dust 
was found to vary by a factor of 12. The column den- 
sity for DR21 (OH ) of 6.3 • 10 17 cm" 2 as determined by 
IZmuidzinas et all is only slightly higher than our result. 

4.1.2. Carbon Monoxide 

The amount of 13 CO is determined in two ways as- 
suming again optically thin emission. Firstly, by using 
T cx ,i3 CO ,„i g h = 79.4K/ln(4/i?i3 CO ,hi g h) (Eq. CUD with 
the 13 CO 6—5 to 3—2 integrated intensity line ratio 
i?i3 COj high an d by Eq. IjOlOp . This ratio is a sensitive tem- 
perature measure from ^20 K up to a few 100 K. The exci- 
tation temperatures vary between 30 and 46 K and repre- 
sent upper limits because of possible self-absorption due to 
a high optical thickness in the lo wer transition. The tem- 
peratu res agree with values from lWilson fc Mauersbergerl 
l|l990|) who find 20 — 65 K from NH3 emission, and are 
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Table 6. Ci and 13 CO line ratio R (in Kkms" 1 ; values taken from Table|3|), excitation temperature (with lower and 
upper lcr limits), and Ci, 13 CO (with relative upper and lower errors), and C 18 column densities (bracketed by T ex 
of 13 CO). The Cn column density was estimated using a method described in Sect. 14. ll All abundances X are relative 
to the H2 column densities (first line; taken from Col. 4 of TableEJ. The last row lists the relative abundance of carbon 
in the three major gas phase species C II, C I, and CO. The CO column density is derived from the lower C 18 column 
density assuming a relative [CO]:[C 18 0] abundance of 490. 



Line/Lineratio 


Source Property 


DR21E 


DR21C 


DR21W 


DR21 (OH) 


DR21 FIR1 


H2,dust 


N(H a ) [10 22 cm" 2 ] 


2.1 


12.1 


3.69 


25.8 


4.84 


Ci 


Rci (Ci 2-1/1-0) 

7ex,[Ci] [K] 

N(Ci) [10 17 cm- 2 ] 
X(Ci) [10" 5 ] 


1.29 ± 0.36 

79lf 

q 9+0.3 

1.5 


1.25 ±0.35 

6-8l° :S 
0.56 


1.06 ±0.30 
56lf 

q Q + 0.3 

1.03 


0.77 ±0.22 

38li 

r 9+O.2 

0.20 


0.99 ±0.28 

ci 76 
01 39 

*-°-0.2 

0.10 


13 CO 6-5/3-2 
13 CO 3-2/1-0 
C 18 1-0 


^3 CO , high ( 13 CO 6-5/3-2) 

r ex ,i3cO,high [K] 

N( 13 CO 6-5) [10 15 cm" 2 ] 
^"co.iow ( 13 CO 3-2/1-0) 

^cx, 13 CO, low [K] 

N( 13 CO 1-0) [10 15 cm" 2 ] 
N(C 18 0) [10 15 cm- 2 ] 
X(C 18 0) [10- 7 ] 


0.62 ±0.26 

43l| 

3.10 ±1.30 

25? 6 7 
33.4i^ 2 9 3 
2.4 - 3.8 
1.2 - 1.8 


0.72 ±0.31 
46g 
jo.o +20 4 

1.59 ±0.67 

1511 
75.6j:^ 2 

8.7-20.0 
0.72 - 1.7 


0.71 ±0.30 
46l 
22.6; 8 5 ;^ 

1.16 ±0.49 

iql6 

40.0t!;J 
3.5 - 8.9 
0.95 - 2.4 


0.27 ±0.12 

soil 

39 6" 10 7 
o».u +21 3 

1.35. ±0.57 
1411 
69.1^: 6 
9.0 - 15.1 
0.35 - 0.59 


0.42 ±0.18 

3530 

97 Q -7 - 2 
z ' - a +13.4 

1.04 ±0.44 

121B 
56.51I; 2 

8.5 - 17.7 
1.7-3.7 


Cn 


N([Cn]) [10 17 cm~ 2 ] 
X(Cn) [10~ 5 ] 


5.25 
2.5 


5.18 
0.43 


3.00 
0.81 


1.47 
0.057 


2.30 
0.48 




[C+]:[C°]:[CO] [%] 


26:16:58 


10:12:78 


16:12:72 


03:10:87 


05:09:85 



always below the C 1 temperatures from above indicat- 
ing an origin from different gas. The 13 CO column den- 
sity, 23 — 54 • 10 15 cm -3 , varies by a factor of 2 as the 
C 1 column density. The same scatter is found if we es- 
timate the column density from 13 CO 1-0 rather than 
13 CO 6—5 (Eq. (|C.11|) ') and derive temperatures by Eq. 
I|C.8|) . Those column densities tend to be higher, between 
33 — 76 • 10 15 cm~ 3 , indicating that a larger amount of gas 
traced by 13 CO is cold. 

Considering the optically thin emission of C 18 0, we 
estimate the total molecular gas column density from Eq. 
HC.12|I using both 13 CO T GX estimates as limits. We find 
high C 18 column densities of about 1 — 2-10 16 cm~ 2 along 
the ridge and 2 — 4 • 10 15 cm -2 toward E and W. 



4.1.3. Ionised Carbon 

In order to get a rough estimate of the C + 
column density from the ISO data, we assumed 
7V(Ch) = 6.4-10 20 xJ( [Cii1) [cm^ 2 ] (ergs" 1 cm- 2 sr" 1 )" 1 
llCrawford et al.lll985h . which is a good approximation 
for an excitation temperature of > 91 K, and a density 
> n cr = 5 • 10 3 cm~ 3 . The peak column density derived 
this way, N(Cn) ~ 5.2 • 10 17 cm~ 3 , is observed towards 
DR21 C and E, while only about 28 % of the peak value is 
found towards DR21 (OH). 



4.1.4. Abundances 

The comparison between C 18 and 13 CO column den- 
sity (for the lower T ox limit in Table 0) shows that the 
abundance ratio [ 13 CO] : [C ls O] = 6.6 - 11.4 is typically 
near the canonical local abundance ratio [490] : [65] = 7.5 
IjLaneer fc PenziasHl990|l . 

The atomic carbon abundances [Ci]/[H2] fall in the 
range of 0.1 — 1.5-10 -5 . The high values are only observed 
at the E and W positions. In the cooler northern positions 
the abundance is reduced by a factor up to 10. Th i s is in 
contrast to the values ofl — 2-10 -5 bv lKeene et al 1 lll997t) 
found in cold dark clouds. 

Except for DR21FIR1, there is a good correlation be- 
tween X(C 1) and X(C 18 0). Toward DR21 (OH), the C 18 
and [C 1] abundance are systematically lower than toward 
DR21 by factors 2 — 3. Although, the position with the 
lowest [Cl] abundance (10~ 6 ), FIR1, has the highest C 18 
abundance (up to 3.7 • 10~ 7 ). 

iMookeriea et al.1 l)2006|) present a compilation of the 
C/CO abundance ratios versus H 2 column densities for 
Galactic star forming regions and diffuse clouds. The ratio 
at the positions DR21E, W, and FIR1, 0.07-0.17, agrees 
well with the relation log[N(C)/N(CO)} = -0MN(H 2 ) + 
19.9 given there for Cepheus B. DR21 and DR21 (OH) 
deviate from this relation but still are within the observed 
scatter. We therefore underline their conclusion that Ci 
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is no straightforward tracer of total H 2 column densities 
and total masses. 

The amount of carbon in the carbon bearing species 
C + , C°, and CO can be compared using their relative 
abundances (last row of TableEJ. The neutral atomic car- 
bon relative abundance is rather homogeneous at about 
9—16 % of the total carbon gas content at all five positions. 
The rest is shared between CO and C II. Ionised carbon 
has the highest abundance toward DR21E (26%). Here, 
the gas is likely exposed to a high UV radiation field and 
carbon species are photodissociated to C I and ionised fur- 
ther to C II. High fractional CO abundances (>85 %), such 
as for DR 21 (OH) and FIR1 wer e also reported for S 106 
C>86 %. ISchneider et all 120031) and W3Main (>60%, 
iKramer et al.ll2004|) . The ratios in other sources are, e.g., 
37:07:56 for IC 63 ll.Iansen et alJll996h . NGC2024 with 
40:10:50 llJaffe fc Plumeill995|) . or r/Carina with 68:15:16 
l|Brooks et al.ll2003|) . 

4.2. Radiative transfer modeling 

In a self-consistent approach, the temperature, column 
and number density, total molecular mass, volume filling 
and mass composition were dete rmined using a 1-D ra - 
diative transfer code (SimLine, lOssenkopf et alJl200l|) . 
In order to describe these physical conditions, we created 
heterogeneous models with internal clumping and turbu- 
lence in a spherically symmetric geometry for each of the 
5 observed regions. The physical parameters of a shell are 
defined by their values at the inner radius and the radial 
dependence is described by a power-law. 

A model with two nested shells provided in most cases 
reasonable results: (i) a hot and dense molecular inter- 
face zone wrapped around a central Hn-region, (ii) en- 
closed by a massive cold cloud - the envelope region. While 
in the inner component primarily the mid- and high-J 
CO emission lines are excited, the envelope is responsi- 
ble for self-absorption and the low lying lines of 13 CO 
and C 18 0. We set the abundance of species with respect 
to molecular hydrogen to standard values (see Table 0) . 
The source size is a crucial parameter that is not well 
confined by the line profiles or the map we obtained. 
A small cloud model can easily produce results similar 
to those for a bigger cloud with a lower volume filling. 
We therefore approximated the r adii using high resolution 
line and continuum obs e rvations tW ils^nfc^laiiexah£rgej| 
1990| lLane et alJll990t Ijaffe et al.lll989t IChandler et alJ 
19931 IVallee fc Fieed l2006j) . The inner cut-off was as- 
sumed to be 0.04 pc (< 10 4 AU), except for DR21 C 
wher e we estimated ~0.12p c from the 14.7 GHz contin- 
uum l|R.oelfsema et al.ll989fl . The Hn- region is considered 
to be free of molecular gas and we neglect the continuum 
contribution by ionised gas. Gradients within the shells 
in any of the parameters are initially set to zero and are 
adjusted during the fit if required. 

We used the following recipe to find appropriate model 
parameters: (i) The linewidth is approximated by adjust- 



ing the turbulent velocity field until both the low-J and 
mid- J line wings match the observed profile. As line emis- 
sion from the inner warm gas is shielded by the colder 
envelope gas closer to the surface, the width of the CO 
absorption dip in the line center constrains the turbulent 
velocity dispersion of this envelope component, (ii) The 
H2-density is defined as a local clump density parame- 
ter and varied between 10 4 and 10 7 cm" 3 independently 
for both shells. The high densities are a pre-requisite for 
high-J CO excitation. (Hi) In the optical thin line wings 
we see predominantly material from the inner shell where 
the simulated macroturbulent motion leads to the higher 
line width which is only partly shielded by the outer shell. 
This allows us to fit the temperatures of both the inner 
and the outer shell from the line shape, (iv) The model 
allows to differentiate between radial and turbulent gas 
motion. Where lines of high opacity toward the line cen- 
ter show a blue-skewed profile, we introduced radial infall 
as a plausible explanation. If possible, we tried to keep the 
radial component close to zero. 

The full line profiles of CO, 13 CO, and C ls O from 
KOSMA and FCRAO are taken into account in a reduced 
X 2 fit in ~1900 degrees of freedom and using 20 — 28 free 
parameters. The second velocity component is masked out 
by a window centered on 9 ± 3 kms -1 . For the ISO lines, 
we fit the integrated flux of three lines. CO 14—13 is not 
considered because no consistent fit was possible in con- 
junction with the other CO lines. In total, we include 11 
lines weighted by the a 2 determined from a baseline re- 
gion not including any lines. The spectra do not contain 
any obvious baseline ripples above the system noise. 

4.2.1. Line Modelling: CO, 13 CO, and C 18 

We independently optimized the model parameters for 
each position. Fig.[3]shows the derived line profiles. At all 
five positions, we obtain solutions of better than \o = 4.1 
which we consider - given the number of fitted lines - a 
very good result, and definitely much better than for a 
single-component model. 

Within the interface-zone, a temperature range from 
80 K (at DR21FIR1) to 148 K (at DR21 C) for gas masses 
of 15 — 60 M Q is found. This temperature is much higher 
than the typical LTE excitation temperature of the mid- J 
CO line peak. The fit converges toward solutions with a 
high density of about 10 6 cm -3 and with a volume filling of 
a few percent, indicating a non-homogeneous distribution. 
The warm gas dominates the emission in the line wings 
and therefore is explained by the steep turbulent velocity 
gradient with velocities up to 30 kms -1 . Because the en- 
velope is much colder and more massive (cf. Sect. I4.2."3|) . 
the average temperature in the regions DR21 C with 36 K 
and DR21 (OH) with 31 K is rather low compared to the 
E, W, and FIR1 position. There, the gas temperature rises 
up to 60 K but the gas masses are only a few hundred so- 
lar masses. Besides the mass fraction, the steepness of the 
gradient in excitation condition and line width between 
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Table 7. The physical parameters for the observed ISO/LWS positions. The radii listed in column 2 give the inner 
and outer edge of each shell. Column 3 to 7 lists the temperature, clump-averaged number density of H 2 , volume 
filling, column density of H 2 , and the gas mass. For these columns, the first entry is the mass-averaged value for the 
whole cloud, values in parentheses show averages within each shell. The last column lists the reduced \ 2 fit result for 
11 lines (6 KOSMA, 2 FCRAO & 3 ISO). 



Position 


Radii 

[PC] 


T gas 

[K] 


n(H 2 ) 
[10 4 cm" 3 ] 


4>v 

volume filling 6 


N(H 2 ) 
[10 21 cm" 2 ]" 


Mass 
[M Q (D/1.7kpc) 2 ] c 


xl 


DR21 E 


0.04/0.2/0.3 


62 (87/55) 


92 (300/4.8) 


0.04 (0.01/0.89) 


46 (9/37) 


243 (46/197) 


2.2 


DR21 C 


0.12/0.2/0.6 


36 (148/32) 


110 (40/110) 


0.024 (0.08/0.02) 


65 (2/63) 


1379 (51/1328) 


3.5 


DR21 W 


0.04/0.1/0.3 


54 (118/52) 


12 (100/9) 


0.59 (0.06/0.8) 


88 (3/85) 


489 (15/474) 


2.8 


DR21 (OH) 


0.04/0.2/0.7 


31 (99/29) 


161 (220/160) 


0.01 (0.01/0.01) 


49 (1/48) 


1421 (34/1387) 


3.9 


DR21 FIR1 


0.04/0.2/0.5 


41 (82/37) 


16 (220/2) 


0.1 (0.01/1) 


47 (4/43) 


700 (59/641) 


4.1 



a Column densities are source averaged and not normalized to 80". 

The product of the volume filling and the abundance of the species is one of the free parameters. We fixed the abundance 

[X/Hs] of 12 CO to 2 ■ 10" 4 (resp. 3.1 ■ 10" 6 for 13 CO, and 4 ■ 10" 7 for C ls O). 
c The molecular gas mass is derived from the cloud-averaged density. This density is the product of local H2 density of the 

clumps and volume filling correction <j> v . 



inner and outer component determines whether we see a 
distinct self-absorption (e.g., DR21 C) or a flat-top (e.g., 
DR21 W) line profile. Table|7|compiles a list of the derived 
physical parameters. A detailed discussion of the individ- 
ual positions is given in App. ^] 

4.2.2. Line Modelling: [Ci] and CS 

As a consistency check, we computed the CS and C I pro- 
files using the same model parameters as derived from CO. 
This implies a homogeneous distribution of CS and C I 
with CO, only affected by the abundances of the species. 

For both species we also obtained excellent fit results. 
The resulting line profiles are shown in Fig. CS 7—6 
is only excited where the H2 density and temperatures 
are sufficiently high. Hence, we find CS emission in the 
dense regions DR21 C and (OH) only. The models fit best 
to the observations f or an abundance of X(CS)=4 • 10~ 9 
l Hatchell et alJ | l998ft. whereas the range 1 • 10" 10 to 2.3 • 
10 -10 given bv lShirlev et al.l (|2003) is clearly too low for 
this region. 

Instead of a carbon abundance variation across the po- 
sitions (as found in Sect. 14.111 . an abundance X(Cl) of 
1 • 10~ 5 nicely fits our [C 1] observations. The C 1 excita- 
tion temperatures as derived in Sect. l4.l"Tl fall in the range 
covered by the two-component radiative transfer model. 

4.2.3. Masses and Volume Filling 

The volume filling, <f) v , is a model parameter capable of 
increasing the clump density without increasing the over- 
all mass as well. The determined <p v is usually low be- 
tween 1% to 8% (Col. 5 in Table GJ. At DR21C, this is 
confirmed by our observations of the CS 7—6 transition 
which is sensitive to densities of several 10 6 cm -3 : Given 
the high H 2 column density (up to 2 • 10 23 cm~ 2 ) this 



translates to a column-to-volume density ratio of about 
0.006 pc. If we assume that the visible extent of the DR21 
source (~0.69pc in CS, after deconvolving the beam) is 
similar to the size along the line-of-sight, a volume fill- 
ing of ^0.01 is derived. A high degree of small scale spa- 
tial and velocity structure is also needed to explain the 
hyper-fine-stru cture anomalie s observed in ammoni a data 
toward DR21 llStutzkil Il98l . IShirlev et all ll2003h . who 
studied CS 5—4 in a sample of high-mass star-forming re- 
gions, derived a volume filling with a median of 0.14. This 
is one order of magnitude higher than our results for the 
dense regions DR21 and DR21 (OH). 

The total mass for all 5 positions is 4.2 • 1O 3 M0, 
and mostly shared between DR21 C and DR21 (OH) with 
~1.4 ■ 10 3 M Q each. In direct comparison to the masses 
derived from dust observation in Sect. I3~3~T1 we find that 
our model gives higher masses (xl.6 — xl.8), except for 
the DR21 (OH) model mass. To compare our mass es- 
timate with masss from the literature, often derived for 
a different distance, we have scaled the mass in Table 
with the square of the distance, basically assuming that 
the column density, computed from the observed lines, 
is independent from the distance. This approach is con- 
sistent with a local radiative transfer treatment, such as 
LTE or escape probability analyses, where the molecu- 
lar excitation does not depend on the total gas mass. For 
the self-consistent radiative transfer code, this assump- 
tion is, however, not completely justified, but deviations 
from this simple scaling relation are to be expected. When 
comparing the values from Table [7| to results obtained 
for much larger distances our quadratically scaled mass 
estimates are therefore always expected to be too high 
and when comparing them to results for lower dista nces 
they should be systematically lower. IShirlev et alJ esti- 
mated virial masses of 561 M Q (D/1.7kpc) for DR21 (resp. 
714M Q (D/1.7kpc) for DR21 (OH)). These estimates are 
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Fig. 8. Cooling curves of the best-fit model line fluxes (marked by open symbols) of CO, 13 CO, and C 18 0. Observed 
line intensities are ind icated as filled symbols (see Table EJ). At position DR21C, we included two observations by 
iBoreiko fc Bet3 l|l99lh of CO and 13 CO J = 9-8. Since the 13 CO 6-5 map does not cover the postitions DR21 (OH) 
and DR21 F1R1 on a fully sampled grid, we expect a larger error. The high- J CO emission at DR21 FIR1 are upper 
limits. 



lower than our masses but were calculated over a much 
smaller radius of 0.17pc(D/1.7kpc ) while our region cov- 
ers up to 0.7 pc. lOssenkopf ei~aT] l|200lh estimated the 
DR21 mass using a radiative transfer analysis on CS data. 
They estimated 578M (D/1.7kpc) 2 at a (cloud-averaged) 
density of 1.6 x 10 3 cm~ 3 . Their lower density is partly 
explained by the larger outer radius, that was chosen to 
be 1.0pc(D/1.7kpc). It is probably also due to the larger 
physical cloud size at the assumed distance of D — 3 kpc 
which prevents the simple quadratic distance scaling of 
the mass in a self-consistent radiative-transfer treatment. 

The total mass of DR21 (OH) corresponds to a cloud- 
averaged column density of ~5 ■ 10 22 cm" 2 . The value of 
5.7- 10 23 cm" 2 reported bv lRichardson et alJ from 1100 
data and our value given in Table [5] based on the far-IR 
dust-emission are higher by a factor of 5 — 10. A flatter 
spectral index (3 may explain why dust and line obser- 
vation ana lyses result in differen t conclus ions. Values be- 
tween 1.0 jVallee fc Fiegell2006|) and 2.0 ()Chandler et all 
Il993h were discussed for DR21 (OH) in the literature. If (3 
is lower than the assumed 1.5, this will have a significant 
effect on the derived A„. For instance, with (3 = 1.25, the 
column density is reduced to about one fourth of our value 
at j3 = 1.5, then is consistent with the CO model results. 

An alternative explanation for the discrepancy could 
be depletion as a resul t of freeze-out of CO onto grains. 
Some authors (e.g., IVallee fc Fieeel I2006L and refer- 
ences therein) use a depletion factor of 10 or more 
(X(CO)=2.5 • 10~ 6 at DR21 (OH)) as app ropriate for star- 
forming clouds. In l.lorgensen et all l(2004l) . the abundance 
is reduced by a factor of 10 for the conditions T < 30 K 
and n > 3 • 10 4 cm" 2 . The DR21 (OH) envelope and the 
resulting column density would comply with these condi- 
tions. However, the CO line data analysis does not allow 
us to draw conclusions toward depletion. 

4.2.4. Problems of the fit 

As noted before, a precise determination of the outer radii 
of the model clouds, and hence their mass, is quite uncer- 
tain and a small increase may already more than double 



the mass. In the models of DR21C and DR21 (OH), we 
found no reasonable solutions for radii much smaller than 
0.6 pc matching the 13 CO and C 18 1—0 observations. As 
possible indication of a systematic error in the column 
density we typically see that the observed peak temper- 
ature of the C 18 1—0 line is weaker than predicted by 
the models (see Fig. while the 13 CO 1—0 line is well 
reproduced. This can either be a result of an inaccurate 
modeling of density and temperature conditions, or our 
assumption about the fixed abundances for CO, 13 CO, or 
C 18 is incorrect. 

The models overestimate the 13 CO 3—2 peak temper- 
ature. Except for DR21 E, this transition is optically thick 
(r > 2) and reflects the kinetic temperature of the enve- 
lope gas. As the CO 3—2 to 7—6 lines also become optically 
thick in the envelope, the brightness temperature at the 
self- absorbing line center of these mid- J lines is governed 
by the same material. The very deep (~5 K) a bsorption 
dip in CO 7—6 observed by Ijaffe et all l)l989|) indicates 
very cold foreground gas (< 20 K) which we cannot re- 
produce with only one component for the outer shell. A 
better match of the observed intensities may be possible 
by adding a third colder layer outside of the envelope. 

Since the cold and warm components certainly are in- 
termixed to some extent rather than separated in two 
nested shells, the focus during the modeling was put on 
a qualitative description of the different components and 
should not be seen as a geometric representation of the 
real cloud structure. We demonstrated, however, that the 
simplified assumption of a spherical cloud is adequate to 
estimate physical properties, such as mass, density, filling 
factors, and temperature. 

4.3. Oxygen and Nitrogen 

4.3.1. ISO/LWS far-IR lines: [Oi], [Om], and [Nil] 

In addition to the high- J CO lines (see Sect. l3~3~2"|) . ISO 

provided the integrated line intensities of [Oi] (63 and 
145 (im) (Table 0J). Their relative strength and cooling ef- 
ficiency depends on the temperature and density regime: 
At low temperatures, rotational mm- and submm-lines of 
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CO dominate the cooling of the molecular and neutral gas. 
At temperatures above T— 150 K, pure rotational transi- 
tions of OH, H 2 0, H 2 , and the 63 /jm [Oi] line become 
more important. [O i] at 63 /jm and the atomic and ionic 
fine-structure lines of [Cn] and [Ci] are diagnostic lines 
for PDRs. Due to the nearby Hn-regions, the ionic lines 
[Om] and [Nil] are also expected to be strong emitters 
and may serve as diagnostics to discriminate the fraction 
of [C n] emission coming from the neutral and ionised gas. 
For completeness, we included the [O in] (52 and 88 /mi) 
and [Nil] (122 /im) lines in the table. 

Since DR21 E, C & W line up along the outflow axis, 
we compare the emission in the outflow lobes with the cen- 
tral position: The strong [O i] emission at 63 /xm decreases 
rapidly toward the west fcf.. lLane et al.lll990tl. This line 
easily becomes optically thick and IPoglitsch et all l|l996[) 
indeed found a deeply absorbed line toward DR21. The 
63 /mi emission extends toward the E position at almost 
the same strength as toward the center. As the A v at the 
same time drops down, the gas traced by Oi must be a 
separate, necessarily warm component. [O i] at 145 /mi is 
peaked at DR21 C, showing only a slight east- west asym- 
metry. 

If we assume thcrmalizcd optically thin [O i] emission, 
we can derive the excitation temperature from T cx = 
227.7K/ln(2. 83/(12.22 •%>!])) (Eq. ftHST l. where R [Ql] 
is the ratio of [Oi] 145/63 /mi emission (cf. Table 0J. The 
ratios of observed intensities, 0.06-0.12, result in temper- 
atures between 170 to 370 K and are to be seen as strict 
upper limits since -ftp i] decreases for optically thick [O i] 
63 /im line emission. The [Oi] line ratio at DR21 (OH) of 
i?[Oi]— 0-29 requires a temperature just below 100 K. We 
detect the weakest emission of [O i] at DR21 FIR1. Errors 
in R[oj] become too large to estimate a temperature here. 

As expected for the position of the Hll-region, the 
ionised [O in] lines both are strong at DR21 C and de- 
creases to the sides. Furthermore, the [O in] 88/52 /mi line 
ratio is a good discriminant of the electron density in the 
range 10 2 — 10 5 cm -3 , while bein g insensitive to a bun- 
dance and temperature variations l|R.ubin et alJll99^) . At 
DR21 C, a ratio of about 0.4 indicates an electron den- 
sity of (N e ) ~ 10 3 cm~ 3 . In the east, the ratio rises to 0.8 
indicating an extension of the H n-region to this side. To 
the west, the ratio is above 1.4. An electron density below 
10 2 cm -3 , consistent with the observed [O in] line ratio at 
DR21 (OH), indicates a low degree of ionisation. No Hn- 
region has been identified in this ar ea, e.g., in the VLA 
6-cm continuum jPalmer et al.ll2004j) . Thus, the finding of 
[O in] emission is a bit surprising. The [O in] r atio > 3±1.3 
found at DR21FIR1 is not covered by the iRubin et alJ 
model but corresponds to a very low (N e ). 

[N n] 122 /im was detected at the E and W posi- 
tions only. At the positon of the Hll-region DR21 C, 
the 3ct noise limit still permits centrally peaked emis- 
sion. The fraction of [C n] stemming from the H n-region 
can be estimated from the intensity ratio of [C n] to 
[Nil]. The Galactic abundance ratio C/N in dense Hn- 
regions leads to an expected ratio of [Cll]/[Nn] = l.l 



((Kramer et alfeooll and references therein) . With the ob- 
served [C ll]/[N n] ratio between 11-19 for DR21, less than 
10 % of the [C n] emission are estimated to come from the 
ionised gas. 

4.3.2. Line Modelling: [Oi] and [Cn] 

A determination of the O I and C II abundances may be 
difficult. If we assume that their abundance is closely con- 
nected to the dense molecular gas, we can calculate and fit 
the emission of these two species by adjusting the abun- 
dance of the two-shell radiative transfer models fSect. l4~2"| . 
At DR21C, we find abundances of 6 • 10~ 5 for Oi and 
2.5 • 10" 4 for Cn. For both Cn and the Oi 63 /im line, 
the opacity toward the line center is about 9 and the line 
profile is deeply self-absorbed. The high energy above the 
ground level of O I leads to almost complete absorption of 
the center emission leaving signatures of warm gas only 
in the wings. The 145 /im transition (r ~ 0.5) is unaf- 
fected and is a better measure of the abundance. As C n 
has only one transition and also shows some absorption in 
the model profiles, the abundance may be less accurate. 
At DR21 (OH), the abundance of O I is -1.2 • 10" 5 and 
9 • 10 -5 for C II. The abundances at all five observed po- 
sitions vary between 1 — 6 • 10~ 5 for O I, and 1 — 4 • 10~ 4 
for C II. 

Compared to our LTE estimate in Sect. 14. ll the Cn 
amount derived here is higher by factors of 13 — 160. 
Because of recombination reactions, it is rather unlikely 
to have a large amount of C n hidden in the cold gas, this 
can be considered as an upper limit. 

5. Line cooling efficiency 

The relative importance of the various line tracers among 
each other and with respect to the FIR continuum pro- 
vides a detailed signature of the heating and cooling pro- 
cesses in the DR21 molecular cloud. Table [S] lists the cool- 
ing intensities at the five selected ISO positions E, C, W, 
OH, and FIR1. 

The total line cooling intensity, /tot, derived by sum- 
ming the contributions of CO, Ci, Cn, and Oi, peaks at 
position C, dropping to less than 15% of this peak value 
at FIR1. 

The relative contributions to /tot vary strongly at the 5 
positions. [Oi] is the strongest coolant at positions E and 
C, while CO is strongest at W, (OH), and FIR1, tying up 
more than 50 % of the total cooling. The relative contri- 
bution of CO to /tot varies between 26 % at E and 76 % 
at (OH) while the relative contribution of [Cn] is lower, 
varying between 7 % at C and 21 % at FIR1. The two [C i] 
fine structure lines contribute generally less than 2 %. The 
ratio of CO cooling relative to [C i] cooling, peaks at C at 
114 and drops to less than 50 at FIR1. 

We also compare /tot with the FIR continuum to derive 
the gas cooling efficiency etot as a measure of the fraction 
of FUV photon energy which goes into gas heating. e tot 
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Table 8. Total cooling intensities of CO, [Ci], [Cn], and [Oi]. For CO we use the best fitting model f Sec. I4.2|) and 
sum over all rotational lines up to J=20. For C I and O I, we sum the observed intensities of both fine structure lines 
respectively. The values in parentheses show the relative contribution to the sum J tot . /jr is the total observed infrared 
intensity (Table [5}. The total cooling efficiency e to t is defined by itot/^iR- For comparison, we also list the e' defined 
by (i[Cn] +/[oi])Mr- 



Intensity in [10 4 ergs 1 cm 2 sr x ] 


DR21E 


DR21C 


DR21W 


DR21 (OH) 


DR21RIR1 


Ico{— E COkj u <2o) 


24.1 (28.8%) 


51.3 (44.4%) 


18.7 (53.7%) 


20.0 (75.9 %) 


13.4 (68.7%) 


/[c.](=£[Ci]) 


0.22 (0.3%) 


0.45 (0.4%) 


0.23 (0.7%) 


0.25 (0.9%) 


0.27 (1.4%) 


7 [0I] (=E [Oi]) 


51.1 (61.1%) 


55.6 (48.2%) 


11.2 (32.1%) 


3.78 (14.4%) 


2.23 (11.4%) 


I{Cn] 


8.2 (9.8%) 


8.1 (7.0%) 


4.7 (13.5%) 


2.3 (8.8%) 


3.6 (18.5%) 


Itot 


83.62 (100%) 


115.45 (100%) 


34.83 (100%) 


26.33 (100%) 


< 19.5 (100%) 


Im 


5500 


33900 


1810 


12700 


3680 


etot [%] 


1.5 


0.34 


1.9 


0.21 


0.52 


^[Cn]+[Oi] ([Cii]+[Oi]63/*m) 


56.4 


57.5 


14.8 


5.24 


5.38 


e' [%] 


1.0 


0.17 


0.82 


0.04 


0.15 



varies by one order of magnitude, between 0.2% at (OH) 
and 1.5 — 2 % at the two outflow positions E and W. 

Table [S] also lists the cooling efficiency e' just taking 
into account the 158 /jm [C II] and the 63 /im [O i] lines, 
ignoring CO, as is often done in the literature, e' varies 
between 0.04% at (OH) and 1% at E. 

Several studies in the literature have determined the 
gas heating efficiency e' in Galactic s tar fo rming regions 
and in external galaxies. IVastel et al 1 ll200lf) find values of 
typically 0.2 — 0.4% in Galactic star forming regions, 1.1 % 
at the Orion Ba r and less than 0.035 % toward W49N. For 
W3Ma.in /TRS5. lKra,mer eTaTI l)2004h rep ort e' of ~ 0.1% 
simila r to the value found in DR21 C. IMalhotra et alJ 
( 200l]) find efficiencies between 0.05 and 0.3 % in 60 unre- 
solved normal galaxies. iKramer et alJ ((2005) find efficien- 
cies between 0.21 and 0.36 % in the centers and spiral arm 
positions of M83 and M51. 

PDR models, e.g., bv iKaufman et all l)l999|) . predict 
the heating efficiency e' of a cloud illuminated from all 
sides between 0.01% to 1.5% and our values are within 
these theoretical limits. 

A remarkable point is that only a few percent of 
the mass (up to 20% at DR21E) is efficiently heated. 
Considering only the cold envelope component of our ra- 
diative transfer models, i.e., removing the warm interface 
zone, shows that the heated gas accounts for most of the 
CO emission (80% to 94%, summing over J up = 1...20) 
in terms of radiated energy. 

Since the total line power of all lines that we consider, 
/tot ) does not include the contributions from other cooling 
lines (e.g., of H2, H2O, OH, and PAHs) we may miss some 
of the cooling energy. 



6. Summary 

We presented the first comprehensive mapping of sub-mm 
lines in the whole DR21/DR21 (OH) region. The line emis- 
sion of [Ci] at 492 GHz and 809 GHz, CO 3-2, 4-3, 6-5, 
and 7—6, 13 CO 3—2 and 6—5, and of CS 7—6 reveals con- 
tributions from different gas components: (i) A very wide 
winged emission in mid-J CO transitions having a line 
width of more than 15 kms^ 1 and self-absorption, and 
(ii) narrow (4 to 6 kms -1 ) emission in the other lines. 
Based on these lines, and supported by mm- and FIR- 
lincs from FCRAO and ISO, we confined temperature, 
density, masses, and volume filling for five positions to- 
ward DR21, DR21 (OH) and DR21 FIR1 using a radiative 
transfer code by modeling the spectral line profiles. 

1. We found that 2-component models consisting of 
a warm inner region and a cold envelope are needed to 
explain the observed strengths of the energetically low 
and high lying transitions of CO and to explain the self- 
absorption. 

2. The bulk of molecular gas at DR21 & DR21 (OH) 
is cold (~30 K) and concentrated in the ridge. The clump 
volume filling is typically a few percent and the density 
within a clump is of the order 10 6 cm~ 3 . 

3. By fitting the dust emission SED, we constrain 
the H 2 column density to 2-5 10 22 cm" 3 for DR21 E, W 
and toward FIR1. At DR21/DR21 (OH) it rises to above 
10 23 cm -3 . The dust temperatures vary between 31-43 K 
and are close to the temperature of the cold gas. 

4. Warm gas (80-150 K) was found at all five positions. 
The warm gas mass is typically a few percent of the to- 
tal mass but accounts for most of the total CO emission 
(80% to 94%, mainly in the mid- and high- J tranitions). 
DR21 (OH), as the coldest source, still shows a significant 
amount of high- J CO emission. 
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5. CS 7—6 emission can be reproduced if the local den- 
sity in the envelope region is sufficiently high. Therefore, 
high densities are expected at DR21 and DR21 (OH) sup- 
porting again the low volume filling factors. 

6. Supported by the model results, we derived the total 
intensity of the strongest emission lines ([Oi], CO, [Cn], 
and [Ci]). Line emission provides up to 2% to the to- 
tal cooling of the atomic and molecular gas and by dust 
grains. At the dense cores DR21 and DR21 (OH), the cool- 
ing efficiency is less than 0.5%. At DR21, Oi and CO are 
equally strong, whereas at DR21 (OH) the CO lines dom- 
inate the cooling with more than 75%. Cn contributes 
generally less than 20% and Ci not more than 1.5% to 
the line cooling. The [O i] emission shows a large variation 
from < 15 % at DR21 (OH) and FIR1 to 61 % at E. 
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Appendix A: Results at individual positions 

Table lists the derived physical parameters (shell- and 
model-averaged) for temperature, H2 number density, vol- 
ume filling factors, H2 column density, and molecular gas 
mass. The resulting spectral profiles are compared with 
the observed lines in Fig. [3] The plots in Figure [H] show 
the fit of the integrated emission for CO, 13 CO, and C 18 
up to J = 20 in comparison to the observed lines. In the 
following, we discuss details of the individual positions: 

DR21C: This position covers the DR21 Hn-region. Its 
surrounding molecular cloud (within a radius of ^0.6 pc) 



forms the southern part of the molecular ridge. Our best- 
fit model indicates that the envelope at a temperature of 
40 K comprises roughly 94 % of the total mass, leaving 6 % 
for the warm (about 150 K) and dynamically active gas 
(wturb betwe en 4.5 and 29 kms -1 ) subje ct to shocks and 
UV-heating. l.Taffe et al.l already found in ll989l indications 
of a warm component with a mass of 55M Q (D/1.7kpc) 2 , 
which is fully consistent with our result from Tabled The 
by far more massive, colder gas component was only indi- 
rectly traced by the CO J— 7— 6 line reversal in their study. 
In our model, the properties of this gas component in the 
envelope shell is much better constrained by the C 18 
and 13 CO J=l— transitions which are close to optically 
thin (r ~ 0.1 resp. 0.8). The turnover at J up pcr ~ 12 (cf. 
Figure |8J results from the increasing critical density lead- 
ing to sub-thermal excitation of the high - J transitions at 
the cl ump density of about 3 ■ 10 6 cm -3 . iBoreiko & Betzl 
l)l99ll) observed CO and 13 CO 9—8 at comparable reso- 
lution onboard the Kuiper Airborne Observatory (KAO). 
We indicate their integrated fluxes in Fig. for DR21 C. 
There is a good agreement with 13 CO, but CO J=9— 8 
is weaker than the here presented J=7— 6 line. This may 
indicate a stronger flattening of the CO rotation curve 
that could only be explained by a discontinuous density 

or temperature structure. 

In support of a shock driven excitation, lLane et "all 
l|l990|) gave an upper limit to CO J=22-21 of 2 • 10" 4 erg 
s _1 cm~ 2 sr -1 in a 44". 8 beam. Considering beam di- 
lution, this limit is still high compared to the inten- 
sity of CO 17—16, the highest rotational line considered 
here. Our prediction with the radiative transfer model of 
2 • 10 _7 erg s _1 cm~ 2 sr _1 can explain why this detection 
attempt failed, albeit we cannot fully exclude the presence 
of > 200 K warm gas. 

DR21 E: The east side of the DR21 Hn-region is not lo- 
cated on the molecular ridge, but enriched with gas ejected 
from the rigde. The lack of cold gas gives raise to an av- 
erage temperature of 62 K for a gas mass of ^240 M Q . 
From the Spitzer/IRAC image it appears that a fraction 
of gas in the east is undergoing a blister phase. Unlike 
the western outflow side, it is not stopped by a dense 
cloud and expands unhampered. We find that 20% of 
the gas is heated to about 90 K and accelerated to ve- 
locities up to v — 32 kms -1 . As the radiative transfer 
code is not designed for outflow geometries, we applied a 
turbulent velocity gradient from 32 to 5 kms -1 to simu- 
late the outflow. For the envelope gas, the asymmetry of 
the CO profiles indicates a motion away from the core 
at ~ 0.7 krns" 1 . Because of the steep drop of high-J 
CO emission, the local clump density was fitted below 
n(H 2 ) = 10 5 cm~ 3 . 

DR21W: Molecular line emission and radiative transfer 
models of DR21W are similar to position DR21E: The 
gas is warmer than in the ridge (> 50 K) but not at densi- 
ties above 10 5 cm -3 . The warm (~ 120 K) gas contributes 
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with a small fraction of about 3.5 % to the total derived 
mass of ~ 430 M . A turbulent velocity gradient between 
30 and 11 kms' 1 is responsible for line wing emission. 
The CO lines show only weak self-absorption although the 
opacity is high at the line center. This is plausible because 
the cloud surface is illuminated by the cluster at DR21 C. 
Hence, we had to implement a positive temperature gra- 
dient from 44 K to 55 K across the cloud envelope. 

DR21 (OH): The position of DR21 (OH) is dominated by 
cold gas from the ridge at a temperature below 30 K over 
a size of up to 0.7 pc. For the modeling, we simplified the 
scenario of sev e ral em bedded cores as reported by, e.g., 
iMangum et alJ (^992) (at possibly different systemic ve- 
locities), to a single clumpy core defined by the warm in- 
terface region with an outer radius of 0.2 pc. The signif- 
icant mid- and high-,/ CO emission is a clear indication 
for warm and dense gas. We find only a few percent (2 %, 
or 34 Mq) of the gas is heated to ^100 K and at densities 
of 10 6 cm -3 or more. There is no evolved o u tflow asso- 
ciated with this positio n, but ICarden et alJ (^986) and 
iRichardson et all lll994h see indications of a very young 
outflow associated with high temperatures and shocks. In 
our model, half of the warm gas mass is at relative ve- 
locities 5<i>turb<32 kms -1 creating the line wing emis- 
sion. The observed blue-skewed CO line profiles are com- 
patible with radial infall motion of the envelope gas of 
1.0-1.5 kms" 1 . 

DR21FIR1: The northern FIR-region contains at least 
three sub-sources. Besides FIR1, the 80" beam partly cov- 
ers the second source FIR2, so we have to consider this as 
contamination. The modeled outer radius overlaps slightly 
with the neighbouring source, but is well within the beam. 
Because high- J CO is not detected at DR21 FIR1, we can 
constrain the temperature to < 82 K. Given the CO 7—6 
flux at this position, the lower limit of the temperature 
in the inner region is T> 50 K. The relatively weak self- 
absorption of the low- J lines is probably due to the pres- 
ence of moderately warm gas (up to 40 K) in the outer 
surface layers, supporting the flat density profile used in 
the model. 

Appendix B: Correction for Reference- Beam 
emission 

Before the implementation of an On-The-Fly-mode 
(OTF), Dual-Beam-Switch (DBS) was the only efficient 
observing mode to use SMART at KOSMA. Therefore, 
observations in 2003 and 2004 of CO 7-6 and [Ci] had 
to be performed in DBS-mode. In this mode, a chopping 
secondary mirror is used for fast periodic switching be- 
tween an ON and an OFF position (Signal- and Reference 
phase). After some time, usually 20 sec, determined by the 
stability of the system, the telescope is moved so that the 
Reference phase is now the ON, and the new OFF position 
falls exactly on the opposite side. In the data reduction, 




SQ^B^O' 00" 00° 



R.A- (J2000) 

Fig. B.l. On-The-Fly map of the 13 CO 2—1 emission 
(SBS2006). An emission- free position is indicated by the 
triangle on the right and was used as OFF position for 
the OTF observations. The box indicates the map size of 
DBS observations in Fig. and The circle in the top left 
corner indicates the 120" beam size. The smaller circles la- 
beled ON1&ON2, OFFl, and OFF2 show an examplary 
DBS measurement (cf. Fig. lB.2(l with a field rotation angle 
a = 10°. 

each OFF position is weighted equally. Slight differences 
in the optical path lengths usually lead to strong standing 
waves in the spectra when using a simple BS-mode. The 
period of these waves corresponds to the optical length to 
the subreflector. These standing waves are very efficiently 
suppressed in DBS-mode. However, mapping of extended 
sources is rather difficult if the chopper-throw (between 
5' — 6') is small compared to the size of the source and the 
chopper movement is fixed to one direction (azimuth for 
KOSMA). Under those conditions, mapping needs care- 
ful planning to prevent self-chopping within the source. 
Observing time is thus not only constrained by source 
visibility but also by sky rotation. For an hour angle of 
h , the secondary throw is only in right ascension, while 
for other angles there is a component in declination. The 
13 CO 3-2 Position-Switched OTF-map in Fig. IB~T1 shows 
that the DR21 ridge is N-S orientated and field rotation 
angles a of less than ±30° are preferred in order to pre- 
vent chopping into the ridge of emission with both OFF 
beams. 

To reduce the impact of artefacts remaining after the 
standard DBS reduction pipeline, all DBS data were cor- 
rected in the following way (see Fig. E21 for an exam- 
ple): (i) Regardless whether the DBS-spectra (Fig. IB. 21 
(a)) show self-chopping or not, the emission in the differ- 
ence (Tig. IB.21 (b)) of the two OFF positions is analysed. 
Any significant signal here is considered as pollution. This 
assumption holds as long as at least one of the two OFF 
positions is free of emission and no absorption against the 
continuum is present. In the example fFig. IB.2l fc) and (d)) 
only the OFF-position of Beam 1 (OFFl) is contaminated. 
(ii) If a signal in the difference is detected above the noise 
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Fig. B.2. This plot demonstrates the method to correct 
self-chopping in Dual-Beam-Switch observing mode. The 
uppermost spectrum (a) shows a standard DBS spec- 
trum of Ci 1—0 with a self-chopped absorption line 
(ON-0.5(OFF1+OFF2)), and (b) shows the difference 
signal (OFF2-OFF1). In the middle, (c) and (d), both 
beams of the DBS spectrum are displayed (Beam 1 or 
ON1-OFF1, and resp. for Beam 2). The result (e) (after 
removing standing waves) shows a faint emission line. 

level, the algorithm tries to reconstruct the emission sig- 
nal (Fig. lB.2l (d)) by removing Beam 1 from the final spec- 
trum (Fig. nOl (r)). i.e. we use only ON2-OFF2 here. This 
scheme works well for situations where standing waves are 
negligible. In cases where standing-wave patterns are per- 
sistent, we transform the spectra to fourier-space, mask- 
out the standing wave frequencies, and transform the data 
back to user-space. Due to the shorter integration time, 
this procedure increases the noise level by a factor of y/2. 
In case of DR21, this method successfully recovered the 
emission because one OFF position was always nearly free 
of emission. If, however, the source extends over more than 
twice the chop-throw in all directions this method does not 
work. 



and 



N 2 = 4.72 10 1 



T mb (Ci 3 P 2 - 3 Pi)dt; . _ 2l 
- cm 



(C.3) 



Kkms 1 

To derive the exitation temperature, we assume LTE 
and substitute N\ and N% using Boltzmann statistics: 

38. 8K 38. 8K 



^-e s= 

9i 



T, 



hv 2 
k B 



(C.4) 
(C.5) 



Rc i is the ratio of integrated intensities (defined in (|C.2|I 
and l|C.3Jl ) and the statistical weight is gj—2J + 1. 
Although the temperature of the 3- level state of [C i] is 
very sensitive to small variations in Rci, the toal column 
density depends only weakly on the temperature: 



2 

N(Ci) = ^Nj 

j=o 
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T mb (Ci 3 Pi- 3 P )dt; . _ 2l 
cm 



K km s 



(C.7) 



where Z = 2_/j=o SJ e fcBTcx is the partition function. 

To extimate the excitation temperature from the two 
transitions 13 CO J=3-2 and J=l-0 we used Eq. (|C.4jl 
which then translates into 

{E 3 -E{)/kB 26.4K 



.low — 



ln(#^) 



In (9/-Ri3 CO ,iow) 



(C.8) 



The transition J— 6— 5 and J=3— 2 gives accordingly 
79.4 K 



IMP V ( C - 9 ) 

In (4/i£i3 CO ,ughJ 

i?i3co,iow (resp. i?i3co,high) is the integrated intensity ra- 
tio of 3-2/1-0 (resp. 6-5/3-2). 

The 13 CO and C 18 column density can then be ap- 
proximated by these equations 

7V( 13 CO) = 1.268- 10 12 (T OX + 0.88 K)e 1111K/T - x 
T mb ( 13 CO(6-5))d« . _ 2 . 



Kkms 1 



(CIO) 



Appendix C: LTE Methods 

A simple determination of the excitation temperature and 
the column density of C I can be done if we assume the 
emission is optically thin and the level popolation is given 
by thermal excitation. The column density of a material 
in an excited state u is given by: 

_ 87rfc B ^ 2 f T mh dv 
u ~ hc 3 A u J Kkms- 1 ' 1 ' 

With Einstein A coefficients for atomic carbon A\ — 7.9 x 



lO^s" 1 and A 2 = 2.7x 10~ V 1 l|Schroder et al.ll99l[) the 
column densities for each transition are 

r m b(Ci 3 P 1 - 3 P )d u _ 2 



Ni = 5.96 10 1 



K kms -1 



(C.2) 



(C.ll) 



A( 13 CO) = 4.71 • 10 13 (T CX + 0.88 K)e 5 - 27K/Tox 
T mb ( 13 CO(l-0))dz; 

Kkms" 1 [ J ' 

7V(C 18 0) = 4.60- 10 13 (T CX + 0.88 K)e 5 27K / T - x 

r mb (c»o(i-o))dt, 2 

Kkms 1 

Using Eq. (|C4|I , an upper limit to the excitation temper- 
ature of Oi can be derived from the [Oi] 145/63 /xm ratio 
(fluxes in ergs -1 cm~ 2 sr" 1 ) by: 
_ 227.7K 
cx " ln(2.83/(12.22-i? [0l] ))- ^ 6) 



